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The influence of uncorrelated (nonmagnetic) overlayers on 
the magnetic properties of thin itinerant-electron fllms is inves- 
tigated within the single-band Hubbard model. The Coulomb 
correlation between the electrons in the ferromagnetic layers 
is treated by using the spectral density approach (SDA). It is 
found that the presence of nonmagnetic layers has a strong 
effect on the magnetic properties of thin fllms. The Curie tem- 
peratures of very thin fllms are modifled by the uncorrelated 
overlayers. The quasiparticle density of states is used to analyze 
the results. In addition, the coupling between the ferromagnetic 
layers and the nonmagnetic layers is discussed in detail. The 
coupling depends on the band occupation of the nonmagnetic 
layers, while it is almost independent of the number of the non- 
magnetic layers. The induced polarization in the nonmagnetic 
layers shows a long-range decreasing oscillatory behavior and it 
depends on the coupling between ferromagnetic and nonmag- 
netic layers. 



I. INTRODUCTION 

In the past few years, there has been growing in- 
terest in the magnetic behavior of thin metalhc films. 
This attention originates both from fundamental physics 
and from applications. The reduced symmetry and the 
lower coordination numbers at the surface offer the pos- 
sibility of inducing new and exotic phenomepa, such as 
perpendicular anisotropy in ultrathin fLlmsEJ,|^interlayer 
couplings and giant magnetoresistance(GMR)Q in ferro- 
magnetic metal/nonmagnetic metal (FM/NM) multilay- 
ers. In experiment, it was shown that ultrathin transi- 
tion metal films can cLisplay long range ferromagnetic order 
from a monolayer onQ. It is well known, that due to the 
Mermin- Wagner theoreniH, an effectively two-dimensional 
spin-isotropic system cannot display long-range magnetic 
, order at any finite temperature. However, in real ultra- 
' thin films, there always exists an anisotropv-. which allows 
almost two-dimensional magnetism to occurcl. 

In experiment, the magnetic thin films arc normally 
grown on a nonmagnetic substrate. The effects of the 
presence of nonmagnetic overlayers on the magnetic prop- 
erties of tiese thin films,— ace. widely studied both in 
experimentCTE^ and theorynlrEj. It has been shown that 
the interfaces between magnetic and nonmagnetic layers 
play an important role with xespect to the magnetic prop- 
erties of multilayers systemsQ. The composition of the in- 
terface between FM and NM has a strong effect on the 
magnetic properties of the FM filma3. The critical Co film 
thickness for the reorientation transition of the magnetiza- 
tion has been observed to shift from 3-4 ML up to 18 ML 



by the presence of carbon at the interfaced. In addition, 
very low coverage of nonmagnetic materials on top of a 
magnetic layers can h^af a strong effect on the direction 
of the magnetizationEiiiS. For example, an almost 90° ro- 
tation from in-plane direction to out-of-plane direction of 
the magnetization is induced by 0.03 monolayer(ML) Gu 
on a 7 ML thick Co film on a stepped Cu(OOl) surfaced. 
The magnetic properties of ultrathin films do not change 
monotonically with the overlayer thickness in some sys- 
tems. The magnetic anisotropy of thin FM film is ob- 
served to rotate from in-plane direction to out-of-plane 
direction with very thin NM overlayers, while it will re- 
turn to-iitplane direction when the NM overlayers become 
thickeiBil3. First principle calculations for a Co ML on 
Cu(lll) show a switch from the in-plane anisotropy of the 
uncovered Co monolayer to perpendicular anisotropy only 
for a 1 ML thick Cu overlayer, pSsihile 2 ML Cu produce a 
slight in-plane anisotropy agaircil. 

The induced polarization in the NM layers has been 
investigated experijaentally in—ipji-iaus systeias- such 
as Pt/f 
Cu/C 



Pd/Niy, Pd/Fctlll^ME3, Pd/CoEHlZl, and 
The direction of the induced polarization 
has been found to be different for different systems. In 
Pd/Fe films, the Pd nearJie interface is ferromagnetically 
coupled with the Fe filnJlj. The induced moment of the 
Pd has been shown to enhance the Curie temperature Tc 
of the filmEJ. On the other hand, in Pt/Co systems with 
thicknesses up to 1.5 ML the Pt layetis negatively spin po- 
larized (antiferromagnetic coupHng)llj. In Cu/Co systems 
the spin polarization of the Cu layers has been shown to 
oscillate as a function of the Cu thickness. This oscilla- 
tion is regarded as the origin of the oscillation Q^i^^ in- 
ter layer-magnetic coupling in Cu/Co multilayeral5 c3. In 
theoryBH, based on first principle calculations, theasplar- 
ization hats been found to be positive both in Pd/F(£3 and 
NM/TMEI (NM denotes noble metal such as Au, Ag and 
Cu, TM: Fe, Cr) systems. Systematic calculations show 
that Fe, Co and Ni overlayers on Pd favour a ferromag- 
netic configuration, whereas V, Cr andJJn overlayers lead 
to an antiferromagnetic superstructuretj. 

Clearly, nonmagnetic overlayers play a very important 
role with respect to the magnetic properties of thin metallic 
films. Theoretically, the influence of nonmagnetic overlay- 
ers on transition metal filpas.is.rnostly investigated within 
first principle calculationscjc^EJ. However, these calcula- 
tions can give the zero temperature properties of the films 
only. Certain idealized model systems have proved to be a 
good starting point for investigating the magnetic behavior 
at finite temperatures. In order to investigate the temper- 
ature dependence of the magnetization, the Ising model 
was adopted in Ref. |2^. However, in transition metals, 
the magnetically active electrons are itinerant. It is by no 
means clear to what extent the results obtained by local- 
ized spin models are applicable to transition metal films. 
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In addition, the quantum interference of Bloch waves in 
the NM layers has been applied to study the-jnterlayer 
magnetic coupling in the FM/NM multilayer£3. Within 
this theory, only the properties of the NM layers are con- 
sidered, while the properties of the FM layers are neglected 
completely. In a theory of interlayer magnetic coupling in 
the FM/NM multilayers, Edwards et al. obtained the dif- 
ference of exchange coupling between two. FM layers for 
difference band occupation in NM layersca. Their results 
corresponds to a Hartree-Fock treatment of a Hubbard- 
like one-band tight-binding model with on-site interaction 
U = in the NM and with f/ = cx) in the FM. However, 
they didn't calculate the magnetization of the FM and NM, 
and the effect of the NM on the properties of the FM was 
neglected. 

The aim of the present paper is to investigate the in- 
fluence of NM overlayers on the magnetism ot|itinerant- 
electron thin films within the Hubbard modelE3. The in- 
fluence of the coupling between the FM and NM layers on 
both the magnetic properties of the FM layers and the in- 
duced magnetization in the NM layers will be investigated 
in detail. 

The Hubbard model was originally introduced to ex- 
plain band magnetism in transition metals and has be- 
come a standard model to study the essential physics 
of strongly correlated electron systems over years, such 
as spontaneous magnetism, metal-insulator transition and 
high-temperature superconductivity. It incorporates in 
the simplest way the interplay of the kinetic energy, the 
Coulomb interaction, the Pauli principle and the lattice 
geometry. In systems with reduced translational symme- 
try, the Hubbard model has been successfully applied to 
describe the temperature-driven reorientation transition of 
magnetization in itinerjarit-electron filmsEJ, metajl-insulator 
transition in thin filmsE2Lthe surface magnetisn£3 and low 
dimensional magnetismEj. The model, though in principle 
rather simple, nevertheless, provokes a non-trivial many- 
body problem, th£ii_could be solved up to now only for 
some special casesETo. In two and three dimensions, 
one still has to resort to approximate treatments. Due 
to the reduced translation symmetry in thin films even 
more complications are introduced. Recently apgencral- 
ization of the spectral density approach (SDA)cjO has 
been applicci t.o j^tiid^t the magnetism of thin metallic films 
and surfacesEElEjO'O The SDA_ which reproduces the ex- 
act results of Harris and Langdlll that concern the general 
shape of the spectral density in the strong-coupling limit 
{U ^ W, U : on-site Coulomb interaction, W: bandwidth 
of the Bloch density of states), leads to rather convincing 
results-cpaoerning the magnetic properties of the Hubbard 
mode]E3c3il3. By comparison with different approxima- 
tion schemes for the Hubbard model as well as numerically 
exact QMC calculaJions in the limit of infinite dimensions 
it has been showiO that the cocrect inclusion of the ex- 
act results of Harris and Langecll in the strong coupling 
regime is vital for a reasonable description of the magnetic 
behavior of the Hubbard model, especially at finite tem- 
peratures. 

The structure of the paper is as follows. First, the 
Hamiltonian of our model is proposed. Then the SDA 



for the Hubbard film is described in a simple way. In sec- 
tion IV we show the results of the numerical evaluation of 
the theory and discuss the magnetic behavior of the film 
system in terms of the layer and temperature dependent 
magnetizations and the quasiparticle density of states. Fi- 
nally, a summary will be given. 



II. HAMILTONIAN OF THE MODEL 

In this paper, we will concentrate on the essence of 
the effect of NM overlayers on the magnetic properties of 
itinerant-electron thin films, as well as the influence of the 
FM layers on the NM layers. The structure discussed here 
is a NM/FM/NM sandwich structure. We study the sym- 
metric situation where the numbers of the NM layers both 
above and below the FM layers are assumed to be equal. 

The description of this film geometry requires some care. 
Each lattice vector of the film is decomposed into two 
parts: 



R, 



(1) 



Ri denotes a lattice vector of the two-dimensional Bravais 
lattice of the surface layer with N sites. To each lattice 
point Ri a c?-atom basis (a = 1, 2, • ■ • , d) is associated, 
referring to the d — 2d]sfM + dpM layers of the film. Here, 
djsfM denotes the thickness of the NM layers and dpM the 
thickness of the FM layers. Within each layer we assume 
translational invariance. Then a Fourier transformation 
with respect to the underlying two-dimensional (surface) 
Bravais lattice can be applied. 

The considered model Hamiltonian reads as follows: 



u 



-Z- ^ V{a)niacrnia- 



(2) 



where cf^^{ciacr) stands for the creation (annihilation) op- 
erator of an electron with spin a at the lattice site Ria, 
niaa = c^o-CiQCT is the number operator, and T^^ denotes 
the hopping-matrix element between the lattice sites Ria 
and Rjf}. The hopping-matrix element between nearest- 
neighbour sites is set to —tpM {—tNu) in the FM (NM) 
layers and to —t^p between FM and NM layers. U denotes 
the on-site Coulomb matrix element and fi is the chemical 
potential. The Coulomb interaction between the electrons 
is only considered in the FM layers. Therefore we choose 
for V(a): 



Via) 



0, a < dMM or a > djMM + dpM, 

1, df^M < a < dNM + dpM- 



(3) 



In the following, all quantities related to the NM layers 
will be labelled by a subscript NM , those related to the 
FM layers by a subscript FM . 
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III. SPECTRAL-DENSITY APPROACH TO THE 
HUBBARD FILM 

Recently a generalization of SDA has been proposed to 
deal with the modifications due to reduced translational 
symmetry. In the following we give only a brief derivation 
of the SDA solution and|i:c|erptliej:eader to previous papers 
for a detailed discussionEJ'LJa^cB. 

The basic quantity to be calculated is the retarded 
single-electron Green function 



(4) 



All relevant information about the system can be obtained 
from the Green function. For example, the spin- and layer- 
dependent quasiparticle density of states (QDOS) is deter- 
mined by the diagonal part of G^'^ (E) : 

The band occupations are given by 

/ + 00 
dE!^{E)p^„{E), 



(5) 



(6) 



where f-{E) is the Fermi function. Fcrromagnetism is 
indicated by a spin-asymmetry in the band occupations 
ria^ leading to non-zero layer magnetizations nia = '^at ~ 
Uai- The band occupation in each layer is given by Ua = 
na-[ + rial- 

The equation of motion for the single-electron Green 
function reads: 



E 



rpO-1 



(7) 



Here we have introduced the electronic self-energy S^*^ (E) 
which incorporates all effects of electron correlations. The 
self-energy will automatically vanish within the NM layers 
due to V{a) = 0. The local approximation for self-energy 
in the FM layers is adopted, which has been tested recently 
for the case of reduced translational symmetryCJ. Because 
the translational invariance is assumed within each layer, 
we then have J:tf,{E) = 6ffi:^{E). 

The key point of the SDA is to find a reasonable ansatz 
for the self-energy in FM layers. Guided by the exactly 
solvable atomic limit of vanishing hopnixig {tpM = 0) 
and by the findings of Harris and LangecZl in the strong- 
coupling limit {U/tpM » 1), a ope-pole ansatz for the 
self-energy can be motivatedE3: 



E 



920 



E 



930 



(8) 



where the spin- and layer-dependent parameters gf„, 92^ 
and 9^^ are fixed by exploiting the equality between two 
alternative but exact representations for the moments of 
the layer-dependent QDOS: 



{m)af3 



1 3 (7 



-in^^lZdEiErGtfAE) 
([[•••[c».,H]----,H]_,c+ J+) 



(9) 



i-times 



Here (• • •) denotes the grand-canonical average and 
[■••,•• •]_^_|_^is the commutator (anticommutator). It has 

been shownH that an inclusion of the first four moments 
of QDOS (to =0-3) is vital for a proper description of fcr- 
romagnetism in the Hubbard model, especially for finite 
temperatures. Furthec- the first four moments represent 
a necessary conditionH to be consistent with the strong 
coupling results of Harris and Lange. Taking into the first 
four moments to fix the three parameters 9 ?^ g^^ . and g^^ 
in (@), one obtains the following self-energyOE3cl: 



E + fl- 



E + p - Ba-o -U{1- Ua^cr) 



(10) 



The self-energy depends on the spin-dependent occupation 
numbers riaa- and the so-called band-shift Bacr that consists 
of higher correlation functions: 



Oaa — -La 



1 



J/3 



Tt/{cta„Cjflai2nia^„ - 1)). 



(11) 



A spin-dependent shift of the band center of gravity 
Ua-aBa-a ((1 — na-a)Ba-a) may generate and stabilize 
ferromagnetic solutions. Although Baa consists of higher 
correlation, fun t. ti pns, it can be expressed exactly via paa 
and E^(£;)BEIM 

1 1 /■+°° 

Baa =Tr + Tfi / dEJ_{E) 

X {^j:-{e -^,)-l)[E- -p)- 7;r] pUe). 

(12) 

Now a closed set of equations is established via Eqs.(||), 
(H), (|l^) and (|l^, which can be solved self-consistently. 



IV. RESULTS AND DISCUSSION 

In our calculations, a fcc(lOO) geometry is assumed for 
both the FM and the NM layers. We consider a uniform 
hopping tpM = tNM = tNF = 0.25eV between nearest 
neighbour sites. The band occupation of the FM layers 
is set to npM = 1-4 for all calculations. The band occu- 
pation in the NM layers is denoted by nj^M- By adjust- 
ing the on-site hopping integrals T^"" we explicitly exclude 
charge transfer within the FM and the NM layers. Further, 
we keep the on-site Coulomb interaction in the FM layers 
fixed at [/ = 12eV, which is three times the bandwidth of 
the three-dimensional fee lattice and clearly refers to the 
strong-coupling regime. In the following we will refer to 
the considered structure as d^M / dpM / dNM- 
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FIG. 1. The layer-dependent magnetization rric, of the FM 
layers for a cLnm / dpM / d-NM sandwich structure as a function 
of temperature for different numbers dpM of FM layers, (a) 
dpM = 1 and 2, (b), (c), (d) refer dpM = 3,4 and 5, respec- 
tively. For (a)-(d), the number of NM layers is d]VA/ = 4. (e) 
dpM = 5 and djvM = (free standing FM film), a = 1 denotes 
the interface FM layer, and a = 2, 3 the inner FM layers. Other 



parameters are: mfm 
Curie temperature. 



i.4, Umm = 0.8, Tg refers to the bulk 



In Fig. 1, the layer-dependent magnetizations are plot- 
ted as a function of temperature for different numbers of 
magnetic layers. The magnetizations are strongly layer- 
dependent. Without NM overlayers (Fig. le; 0/5/0 struc- 
ture), the magnetization in each layer is fully polarized 
at very low temperatures. The magnetization curves of 
the inner layers show the usual Brillouin-type behavior. 
The surface magnetization, however, shows a very differ- 
ent behavior: It depends almost linearly on temperature 
in the temperature range T/T(f = 0.7-0.9. Compared to 
the inner layers, the surface magnetization decreases sig- 
nificantly faster as a function of temperature and shows a 
tendency to a reduced Curie temperature. However, due 
to the coupling between the surface and the inner layers 
that is induced by the electron hopping, aj-unique Curie 
temperature for the whole film is obtainedEj. When the 
NM overlayers are taken into account (see Fig. la, lb, Ic, Id; 
4/ciFM/4 structure), the interface magnetization of the FM 
layers and its temperature behavior are strongly affected 
by the interaction between the FM and NM layers which is 
induced by the electron hopping at the interface. The in- 
terface is no longer fully polarized at low temperatures and 
decreases more roundly as a function of temperature than 
in the case without NM overlayers. The linear dependence 



of the interface magnetization in the range of T/T^ = 0.7- 
0.9 has disappeared. The magnetization of the inner FM 
layers is only weakly affected by the presence of the NM 
layers (see Fig. Id). 
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FIG. 2. The Curie temperature Tc as a function of the num- 
ber dpM of FM layers for different numbers djvA/ of NM layers. 
djvM = refers to the free standing FM film. The parameters 
are: npM = 1.4, njvM = 0.8. 

From Fig. 1, one can also read off that the Curie tem- 
perature Tc increases gradually as a function of the num- 
ber of the FM layers clpM and reaches its bulk value for 
dpM = 8-10. The influence of the NM layers on the Curie 
temperature is analyzed in more detail in Fig. 2 where 
Tcidpivi) is shown for different numbers of the NM layers. 
As can be seen in Fig. 2, the behavior of Tc as a function 
of dpM is quite different for d^M = and d^M 7^ 0. In 
the case without NM overlayers {d^M — 0) Tc increases 
very steeply a^-iS function of dpi\i and saturates already 
for dpM =3-5E3. The influence of the NM layers is most 
important for very thin magnetic films dpM < 6. For 
dpM — 3,4 the presence of only one NM overlayer leads 
to a reduction of the Curie temperature of about 5% com- 
pared to its free FM film value. It is interesting to note 
that there are only minor changes in Tc when the number 
of NM layers is further increased. Note that for dpM > 3, 
the reduction of the Tc is strongest for just one NM layer. 
As df^M is further increased, the reduction of the Tc be- 
comes slightly decreasing again. Tc as a function of difM 
saturates already for d^M = 3. This indicates that the 
influence of a very thin NM topping on the magnetic prop- 
erties of the FM filmS|-is-S|tronger than for thick NM over- 
layers. In experimentLfO, it has also been observed that 
very thin NM layers have a stronger effect on the magnetic 
properties than thick overlayers. For example, the direc- 
tion of the magnetization of a Co film shows a rotation 
from the in-plane direction to the out-of-plane direction 
induced only by very thin coverage of Cu, while for thick 
coverage of Cu, the direction of magnetization of the FM 
films will turn back to in-planc directionBila. First princi- 
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pie calculations for a Co ML on Cu(lll) predict a switch 
from the in-plane anisotropy of the uncovered Co mono- 
layer to perpendicular anisotropy only for a 1 ML thick 
Cu overlayer, w]fiije 2 ML Cu produces a slight in-plane 
anisotropy agair 




FIG. 3. The Curie temperature Tc and magnetizations iric, 
of a (Inm /dpM /dNM = i/'i/l sandwich structure as functions of 
the band occupation umm of the NM layers, (a) Tc, (b) the 
magnetization of the NM layer, and (c) the magnetizations of 
the FM layers, a — 1 indicates the interface FM layer. The fully 
polarized magnetizations of the FM layers are 0.6. Tq denotes 
the Curie temperature of the free FM film. Other parameters 
are: ufm ~ 1-4. For (b) and (c) the temperature is T = O.ITq. 

In order to investigate the influence of the NM overlayers 
on the magnetism of the FM films in more detail, we have 
calculated the Curie temperature as a function of the band 
occupation unm of the NM layer for a 1/3/1 snadwich 
structure (Fig. 3a). Tc does not change monotonically 
as a function of njvA/- For very small hnm it decreases 
rapidly from the free FM film value T^ {cInm = 0) to a 
minimum at tinm =0.35-0.4. Then it increases from the 
minimum to a maximum (at njvM ~ 1.97) which lies above 
Tq. Finally, it drops quickly to T^ at tlnm = 2.0. The two 
limiting cases (Tc = TjS at n^M = and njvM = 2.0) are 
easy to understand. In these two situations, the NM bands 
are either empty or fully occupied and, therefore, have no 
influence on the properties of the FM film. To understand 
the behavior of Tc as a function of n^M, we have also 
calculated the magnetizations of the NM layer and the FM 
layers as a function of njvA/ at low temperature in Fig. 3b 
and Fig. 3c. The magnetization of the NM layer shows 
a similar dependence on unm as the Curie temperature 



(see Fig. 3b). Since the magnetization of the FM layers 
is always positive, the sign of magnetization of the NM 
layer represents the coupling between FM and NM at the 
interface. A positive sign corresponds to ferromagnetic 
coupling, a negative sign to antiferromagnetic coupling. 
The results indicate that Tc is strongly affected by the 
coupling of FM and NM. From Fig. 3c, one can see that 
the coupling of FM and NM has also a strong effect on the 
magnetizations in the FM layers, especially on the interface 
magnetization at low temperature. 
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FIG. 4. Layer-dependent QDOS of a 1/3/1 sandwich struc- 



ture for difi^erent unm at low temperature (T 
nNM ~ 0.4, (b) riNM = 0.8 and (c) unm = 1.4. 
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Further insight about the coupling between the FM and 
the NM layers can be obtained from the layer-dependent 
QDOS. In Fig. 4 the QDOS of a 1/3/1 sandwich structure 
is shown for T = O.IT^ and three different values of the 
band occupation of the NM layers [tinm = 0.4, 0.8, 1.4). 
For Unm — 0.8 the temperature dependence of the QDOS 
is plotted in Fig. 5 for T = 0.1, 0.9, l.OTc. Two kinds 
of splittings are observed in the FM spectrum. The strong 
Coulomb interaction between the electrons leads to a split- 
ting of the spectrum into a high and a low energy part 
("Hubbard splitting"). These two quasiparticle subbands 
( "Hubbard bands" ) are separated by an energy amount of 
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the order U. In the lower subband the electron mainly 
hops over empty lattice sites, whereas in the upper sub- 
band it hops over lattice sites that are already occupied 
by another electron with opposite spin. The weights of 
the subbands scale with the probability that a propagat- 
ing electron meets the one or the other situation. The total 
weight of the QDOS of each layer is normalized to 1. In 
the strong coupling limit the weights of the lower and the 
upper subbands are given by (1 — Ua-a) and ria-cr, respec- 
tively. Due to the vanishing Coulomb repulsion the Hub- 
bard splitting disappears in the NM spectrum. For temper- 
atures below Tc, an additional spin splitting ("exchange 
splitting") in majority- (ct =t) and minority- (cr =|) spin 
direction occurs in both the FM and the NM spectrum, 
leading to a non-zero magnetization ma = Uai —n-ai. Note 
that, in principle, the NM and the FM spectrum for each 
spin direction occupy exactly the same energy region, thus 
preventing the electrons to be trapped. However, the cor- 
responding spectral weight may become very small as can 
be seen in Fig. 4 and Fig. 5. 

First we want to discuss the riArM-dependence of the 
QDOS at low temperatures. For the case of dj^M = 0, the 
majority QDOS lies completely below the chemical potep^ 
tial, and the system is fully polarized {na\ = 1) at T = OEj 
(see also Fig. le). If the NM overlayer is taken into ac- 
count (Fig. 4), the majority QDOS of both the interface 
and inner FM layers have tails above the chemical poten- 
tial due to the hybridization between FM and NM layers. 
As a consequence the system is no longer fully polarized at 
low temperatures (see also Fig. la, lb, Ic, Id). The weight 
of the tail determines the reduction of the magnetization 
compared to the fully polarized state. Since the weight 
of the tail becomes smaller as the band occupation of the 
NM layer is increased from n^M = 0.4 to unm = 1.4 the 
reduction of the magnetization in the interface gets weaker 
as well (see also Fig. 3c). 

Let us now turn to the NM layer. As mentioned above, 
the QDOS of the NM layer has just one kind of splitting 
- the exchange splitting. The QDOS of the NM layer is 
quite different for majority- and minority- spin direction. 
Because there is no Coulomb interaction in the NM layers, 
the majority- and minority-QDOS do not affect each other. 
It is interesting to note, that above the chemical potential 
the majority NM-QDOS resembles the BDOS (Bloch den- 
sity of states) of the square lattice which is equivalent to 
the free standing fcc(lOO) monolayer. This is because for 
energies above the chemical potential a (cr =t)-electron 
within the NM layer is effectively isolated since the spec- 
tral weight of the (cr =t)-FM-QDOS is very small in this 
energy region. 

For Unm — 0.4 (Fig. 4a), there is very low spectral 
weight in the majority NM-QDOS below the chemical po- 
tential. As a consequence the number of spin-up electrons 
is smaller than the number of spin-down electrons in the 
NM layer. The magnetization of the NM layer is neg- 
ative, i.e., the NM and FM layers are antiferromagneti- 
cally coupled. With increasing unm the center of gravity 
of the QDOS of the NM layer gets shifted to lower ener- 
gies. When the peak of the majority NM-QDOS crosses 
the chemical potential, the number of the majority-spin 



electrons increases faster as a function of Unm than the 
number of minority-spin electrons. The magnetization of 
the NM layer increases and becomes positive (see Fig. 4b 
and 4c, also Fig. 3b). Of course, as unm increases, the 
shape of the QDOS will also change. Until njvM — 1.97, 
the magnetization of the NM layer will increase as a func- 
tion of riNM- Then it drops quickly, because both the 
majority- and minority-spin QDOS gets shifted below the 
chemical potential. 
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FIG. 5. Temperature dependence of the layer- dependent 
QDOS for a 1/3/1 sandwich structure, (a) T = O.IT^, (b) 
T = 0.9Tc and (c) T = Tc. 

The temperature dependence of the FM-QDOS (Fig. 5) 
is dominated by two correlation effects. As a function of 
increasing temperature the spin-splitting between the cen- 
ters of gravity of the majority and minority quasiparticle 
subbands decreases. In addition there is a temperature- 
dependent transfer of spectral weight between the lower 
and the upper quasiparticle subband. The interplay of 
these two correlation effects leads to a rather rapid demag- 
netization of the system as a function of temperature and 
allows for Curie te pijpe |ijatures in a physically reasonable 
order of magnitudeE3c3a. At T = Tc^ the spin-splitting 
has disappeared completely, both in the FM and the NM 
spectrum (Fig. 5c). 



6 



Due to the coupling between FM and NM layers, 
the NM-QDOS becomes temperature dependent as well. 
While the band edges of the NM-QDOS stay fixed, there 
is a redistribution of spectral weight within the NM-QDOS 
as a function of temperature. We would like to point out 
that this temperature-dependence may even result in a 
temperature induced change from ferromagnetic to anti- 
ferromagnetic coupling between FM and NM layers. From 
Fig. 5b, one can read off that at T = O.QTc the number of 
spin-up electrons is smaller than the number of spin-down 
electrons. Due to the reduced spin-splitting in the FM 
spectrum, the main peak of the majority NM-QDOS gets 
shifted above the chemical potential. As a consequence 
the coupling between the FM and the NM layer at the 
interface changes from ferromagnetic (Fig. 5a) to antifer- 
romagnetic (Fig. 5b) coupling. This kind of temperature 
induced change of the coupling is observed only for very 
thin NM layers {dNM = 1). For (Inm > 1, we have not 
found such a behavior. 
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FIG. 6. The polarization of the NM layers for different num- 
bers dNM of the NM layers, (a) Magnetizations of the NM 
layers as a function of the layer index, (b) Magnetization pro- 
file of the NM layers for different band occupations of the NM 
layers, (c) Magnetizations of the NM interface and the surface 
layer as a function of cLnm- Other parameters are: dpM = 4, 
T = Q.\TS. 

Finally, we want to discuss the induced polarization in 
the NM layers which is shown in Fig. 6. The induced mag- 
netization of ±he NM layers oscillates with the layer index 
for cLnm > 4.e3 As a consequence, the total magnetization 
of the NM layers oscillates as a function of the thickness 



of the NM as well. An oscillatory behavior of total magne- 
tization of the NM layers as a function (if-|the thickness of 
the NM has also been obtained by BruncO within the the- 
ory of intcrlayer magnetic coupling. However, our starting 
point and our results are completely different. The quan- 
tum interference was introduced to understand the prop- 
erties of the NM layers. The Bloch waves of the NM layers 
are regarded as to be confined within two FM perturba- 
tions (or one FM perturbation and one Vacuum). Due to 
the constructive and destructive reflection at this confine- 
ment, the densities of spin-up and/or spin-down electrons 
are found to be oscillating as a function of the thickness of 
the NM. While the long range magnetic order was excluded 
in his theory. Within this concept, one can find that the 
magnetization of the NM interface should oscillate around 
zero. On the contrary, in our calculation, the amplitude of 
the oscillation decreases from the interface to the surface 
(Fig. 6a) and the profile of the induced magnetization is 
almost independent of djvAf(Fig. 6a). The induced mag- 
netization at the interface exhibits only a very small os- 
cillation around a certain negative value as (Inm > 4 (see 
Fig. 6c). This indicates that the coupling between FM and 
NM layers is only affected by the properties of the FM and 
NM layers close to the interface, whereas it is independent 
of the thickness of the NM overlayers as cInm > 4. 

Due to the oscillatory behavior of the induced magne- 
tization, wc find the surface magnetization to be also os- 
cillating as a function of the number of the NM layers 
(Fig. 6c). This surface magnetization may be observed 
in experiment by use of spin polarized photo-emission for 
very short mean free paths. We would like to point out 
that the oscillation of the NM magnetizations with respect 
to the layer index is affected by the band occupation of 
the NM layers. The amplitude as well as the period of the 
oscillation changes when the band occupation of the NM 
layers is changed (see Fig. 6b). 



V. CONCLUSIONS 

The effect of uncorrelated layers on the magnetic proper- 
ties in thin Hubbard films is studied by the spectral density 
approach (SDA). The SDA, which reproduces the exact 
results of i/C/-perturbation theory of Harris and LangecJ 
in strong-coupling limit, leads to rather convincing re- 
sults cojaeaHiipff the magnetic properties of the Hubbard 
modelE^HH^L^. By comparison with different approxima- 
tion schemes for the Hubbard model as well as numerically 
exact QMC calcidation in the limit of infinite dimensions it 
has been showrc3 that the correct inclusion of the exact-re- 
sults of Harris and Lange in the strong coupling regimdlJ is 
vital for a reasonable description of the magnetic behavior 
of the Hubbard model, especially at finite temperatures. 

Within our theory, FM and NM layers are treated on the 
same footing. The effects of the coupling between FM and 
NM on the magnetic properties of the FM layers, as well as 
on the NM layers, have been studied in detail. The Curie 
temperature of thin FM films is modified by the presence 
of NM layers. For npM = 1-4 and njvM = 0.8 the Curie 
temperature increases gradually as a function of the num- 
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ber of FM layers and converges to the corresponding bulk 
value for dpM = 8 — 10. The induced polarization in the 
NM layers displays a long-range decreasing oscillation with 
respect to the layer index of the NM layers. The induced 
magnetization of the interface NM layer hardly changes 
as a function of number of the NM layers when d^M > 4. 
This means the coupling between the FM and NM layers is 
determined by the intrinsic properties of the FM and NM 
layers, such as the band occupation, the in-site Coulomb 
interaction and the hopping, and it-is independent of the 
numbers of the FM and NM layersEll. 

The magnetic properties of this thin film system have 
been microscopically understood by means of spin- layer- 
and temperature- dependent quasiparticle density of states 
(QDOS) for a 1/3/1 sandwich structure. There exist two 
correlation induced band splittings in the FM spectrum. 
Besides the Hubbard splitting there is an additional ex- 
change splitting for temperatures below Tc- Due to the 
vanishing Coulomb repulsion the Hubbard splitting disap- 
pears in the NM spectrum. For the case of Unm — 0, 
the majority FM-QDOS lies completely below the chemi- 
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cal potential, and tip*, system is fully polarized {uai 
at zero temperatureEj. If the NM layers is taken into ac- 
count, the majority QDOS of both the interface and inner 
FM layers have tails above the chemical potential due to 
the hybridization between the FM and NM layers. As a 
consequence the system is no longer fully polarized at low 
temperatures. There is a reduction of the FM magnetiza- 
tion compared to the fully polarized state. The reduction 
of the magnetization of the interface FM layer gets weaker 
as the band occupation of the NM layer is increased from 
nNM = 0.4 to riNM = 1-4. In addition for unm = 0.4, 
there is very low spectral weight in the majority NM- 
QDOS below the chemical potential. As a consequence the 
number of spin-up electrons is smaller than the number of 
spin-down electrons in the NM layer. The magnetization 
of the NM layer is negative, i.e., the NM and FM layers are 
antiferromagnetic coupled. As n^M increases, the magne- 
tization of the NM layer increases and becomes positive. 
The change of the coupling between the FM and NM layers 
can also be induced by the temperature. This temperature 
induced change of the coupling is observed only for very 
thin NM layers. 
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